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Abstract 

We study the tunneling properties of a cigar-shaped Bose-Einstein conden- 
sate by using an effective ID nonpolynomial nonlinear Schrodinger equation 
(NPSE). First we investigate a mechanism to generate periodic pulses of co- 
herent matter by means of a Bose condensate confined in a potential well 
with an oscillating height of the energy barrier. We show that is possible to 
control the periodic emission of matter waves and the tunneling fraction of 
the Bose condensate. We find that the number of emitted particles strongly 
increases if the period of oscillation of the height of the energy barrier is in 
parametric resonance with the period of oscillation of the center of mass of 
the condensate inside the potential well. Then we use NPSE to analyze the 
periodic tunneling of a Bose-Einstein condensate in a double-well potential 
which has an oscillating energy barrier. We show that the dynamics of the 
Bose condensate critically depends on the frequency of the oscillating energy 
barrier. The macroscopic quantum self-trapping (MQST) of the condensate 
can be suppressed under the condition of parametric resonance between the 



frequency of the energy barrier and the frequency of oscillation through the 
barrier of the very small fraction of particles which remain untrapped during 
MQST. 
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I. INTRODUCTION 



In the last few years macroscopic quantum tunneling (MQT) with dilute Bose-Einstein 
condensates of alkali-metal atoms has been the subject of many theoretical [1-4] and ex- 
perimental [5-7] studies. In a recent paper we have suggested a mechanism to generate 
pulses of coherent matter by using MQT of a falling and bouncing condensate [8]. In the 
first part of this paper we propose another, simpler, mechanism to produce periodic atomic 
wave-packets by means of MQT of a Bose condensate through the barrier of a potential 
well. In our theoretical investigation we consider cigar-shaped Bose condensates so that 
we can use an effective ID nonpolynomial nonlinear Schrodinger equation, which we have 
shown to be quite accurate when the aspect ratio is larger than 3 [9,10]. We consider a 
Bose-Einstein condensate confined in the vertical axial direction by two Gaussian optical 
barriers and in the transverse direction by a magnetic or optical harmonic potential. By 
periodically changing the height of the lower-lying Gaussian barrier it is possible to generate 
periodic waves of coherent matter: only when the height of the energy barrier is sufficiently 
small a consistent fraction of condensed atoms can tunnel. We show that it is possible to 
control the period of emission of the matter waves and the tunneling probability. We find 
that the emission probability can be strongly enhanced by using the parametric resonance 
between the oscillation of the height of the energy barrier and the oscillation of the center 
of mass of the condensate inside the potential well (resonance- induced coherent emission). 

In the second part of the paper we study another problem of tunnelling, the one of trapped 
particles in a double potential well. In particular, we investigate the effect of a periodically 
varying barrier in the tunneling of the Bose condensate. The periodic oscillations of a 
Bose-Einstein condensate in a double-well has been the subject of many papers [1-4] but, 
to our knowledge, the role of a time-dependent energy barrier has not been analyzed. We 
investigate the problem using our effective ID nonpolynomial Schrodinger equation and also 
the two-mode classical-like equations introduced by Smerzi et al [1]. The numerical results 
obtained with the two methods are in quantitative agreement. In particular, we find that 
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the phenomenon of macroscopic quantum self-trapping of the condensate can be controlled 
(reduced or suppressed) by the frequency of oscillation of the double-well energy barrier 
again via parametric resonance. 



II. MAPPING FROM 3D TO ID 

It is well know that the 3D Gross-Pitaevskii equation (GPE) is accurate in describing 
the Bose-Einstein condensate of a dilute boson gas at temperatures well below the transi- 
tion temperature, where thermal excitations can be neglected [11]. Often the condensates 
produced in experiments are cigar-shaped and this suggests to map the 3D GPE into an 
effective ID equation, which simplifies greatly the solution of the equation. This prob- 
lem is not trivial due to the nonlinearity of the GPE. By using a variational approach, we 
have recently derived an effective time- dependent ID nonpolynomial nonlinear Schrodinger 
equation (NPSE) which describes very accurately cigar-shaped Bose condensates [9,10]. We 
assume that the condensate is confined by a magnetic (or optical) harmonic potential with 
frequency cu± and harmonic length a± = {Ti/muj^) 1 / 2 in the transverse direction and by 
a generic potential V(z) in the axial one. The total wave function of the condensate is 
ijj(x,y, z,t) = f(z,t)g(x,y,t), where the transverse wave-function g(x,y,t) is a Gaussian 



with the width i] such that rf = a\^j 1 + 2a s N\f\ 2 with a s the s-wave scattering length [10]. 
The NPSE is given by 



<4/ 

dt 



h 2 d 2 Jr , . gN 
+ V(z) + 



2m dz 2 2na 2 ± ^Jl + 2a s N\f\ 2 

/, (1) 



2 \^/l + 2a s N\f\ 2 V 

where f(z, t) is the axial macroscopic wave function of the condensate, g = AnTi 2 a s /m is the 
scattering amplitude and iV is the number of condensed bosons and the function f(z,t) is 
normalized to one. 

We observe that in the weakly-interacting limit a s N\f\ 2 « 1, NPSE reduces to a ID 



GPE with the nonlinear coefficient g' given by g' = g/(2na\). Instead, in the strongly- 
interacting limit, NPSE becomes a nonlinear Schodinger equation with the nonlinear term 
proportional to \ f\f- 

In paper [10] we have tested the accuracy of the NPSE in the determination of the ground- 
state and collective oscillations of the condensate with axial harmonic confinement and also 
in the description of tunneling through a Gaussian barrier. In particular, we have compared 
NPSE with the the full 3D GPE. The conclusion is that NPSE is very accurate when the 
aspect ratio in larger than 3 and gives better results than all other effective approaches 
recently proposed [12,13]. 

III. PERIODIC EMISSION OF MATTER WAVES 

In a previous paper [8] we have analyzed the pulsed quantum tunneling of a Bose con- 
densate falling under gravity and scattering on a Gaussian barrier that model a mirror of 
far-detuned sheet of light. In this section we study another mechanism which produces pe- 
riodic emission of matter waves by means of the quantum tunneling of a Bose condensate. 
We consider a cigar-shaped condensate confined by a harmonic potential in the transverse 
direction and under the action of the following external potential in the vertical axial direc- 
tion: 

V(z) = V 1 e -(^i) 2 /- 2 + V 2 e -(^) 2 /- 2 + mgz . (2) 

Thus, in addition to the gravity potential mgz, there are two Gaussian functions that model 
a confining potential well. The shape of this potential is shown in Figure 1. Such a configu- 
ration can be experimentally obtained by using two blue-detuned laser beams (perpendicular 
to the axial direction) which are modelled by the two Gaussian potentials. By varying the 
intensity of the lower-lying laser beam one controls the height of the energy barrier and 
therefore the tunneling probability. In particular we consider the case in which the height 
V2 of the lower barrier is a sinusoidal function of time. As an example we consider a Bose- 
Einstein condensate of 23 Na atoms with a scattering length a s = 30 A. For the frequency oj±_ 
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of the axial harmonic confinement we choose u± = 2rr kHz. Setting oj z = u±/10, we measure 
lengths in units a z = (h/muj z ) l l 2 , time in units uj~ 1 and energy in units Tiw z . In these units 
the period of small oscillations around the minimum of the potential well (2) close to one. 
The NPSE is integrated by using a finite-difference predictor-corrector method. The initial 
wave function of the condensate is found by solving the equation with imaginary time. 

First we calculate the "ground-state" of the system by choosing V\ = 1000, Vi = 900, 
z\ = 17, z 2 = 24 and a = 1.6. Then we reduce the energy barrier V\ studying the time 
evolution and tunneling of the condensate through the potential well. The results are shown 
in Figure 2. There is a periodic tunneling of matter waves and the period of tunneling is 
proportional to the period of small oscillations of the center of mass of the condensate around 
the minimum of the axial external potential V(z). To compare the chemical potential of 
the Bose condensate with the energy barrier Vi, the effective chemical potential /i e ff of the 
Bose condensate is obtained from the full one by subtracting not only the transverse energy 
per particle but also the gravitational energy per particle. 

A simpler mechanism to generate periodic pulses is to have the lower-lying Gaussian 
barrier with a periodic height given by: 

2ir 

V 1 (t) = V a + V b sm(—t), (3) 

T 

where r is the period of oscillation. We choose V a = 300 and V b = 200. In Figure 3 we plot 
the axial density profile p{z) of the matter-waves coming out from the potential well with 
r — 1.5. Figure 3 shows that, as expected, the emission of pulses is periodic. 

We find that both the period of emission and the tunneling fraction Pt depend on r. 
For very small values of r, the fraction of ejected particles is quite small because the height 
of the energy barrier oscillates very quickly. On the other hand, for large values of r the 
fraction P T is again very small. In the inset of Figure 4 the tunneling probability Pt is 
shown as a function of r for a Bose condensate with iV = 10 3 and iV = 10 4 atoms. The 
emission probability has its absolute maximum near r = 0.8, that is the average period of 
the oscillations of the center of mass of the condensate inside the potential well. 
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The emission probability Pt grows by increasing the number N of atoms and it is strongly 
enhanced by setting the oscillating barrier at the resonance condition. Nevertheless, by 
plotting the ratio Pt/Pt as a function of r, where P T is the tunneling probability with 
Vi(t) = V a — Vb = 100, one finds that the resonance is suppressed by the interaction, i.e. by 
increasing N (see Figure 4). Note that in the non-interacting case (N — and /j, e ff = 21.88) 
one finds P° = 1.11 x 10~ 4 , but Pt = 0.35, thus the tunneling ratio Pt/Pt is larger than 
3 x 10 3 . 

The enhancement of the emission probability of coherent matter near r = 0.8 has a 
classical explanation. It corresponds to the condition of resonance between the oscillation 
of the energy barrier and the oscillation of the center of mass of the condensate. Such an 
effect is clearly shown in Figure 5, where we plot the axial coordinate and the phase-space 
portrait of a classical particle under the action of the potential (2). The dynamics of the 
classical particle is obtained by solving the Hamilton equations with a 4th-order Runge- 
Kutta method. In Figure 5, for r = 0.75 the system is close to the resonance condition: 
the periodic energy of the barrier is pumped into the energy of the classical particle and 
the particle is able to escape from the trap. Far from the resonance condition, the periodic 
energy of the barrier is not fully transferred to the energy of the classical particle, so the 
particle remains confined. Close to the resonance condition the phenomenon of coherent 
emission of atoms is no more due to quantum tunneling but to the parametric resonance 
[14] between the period of oscillation of the center of mass of the condensate and the period 
of oscillation of the energy barrier. 

Because of the presence of an oscillating barrier of period r, the period f osc of oscillation 
of the center of mass of the condensate follows the law f osc = r osc (l + e sin(27r£/r)), where r osc 
is the period of oscillation of the center of mass in the absence of an oscillating barrier and 
e represents the amplitude of the periodic perturbation. Text-books [14] predict parametric 
resonance for r = (n/2)r osc with n an integer number. In fact, as shown in Figure 4, 
there are other local maxima in the plot of the emission fraction whose positions follow this 
resonant condition formula. 
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IV. PERIODIC TUNNELING IN A OSCILLATING DOUBLE- WELL POTENTIAL 



In this section we investigate the dynamics of a Bose-Einstein condensate of 23 Na atoms 
in a double-well trap given by a harmonic anisotropic potential plus a Gaussian barrier along 
the z axis, which could model the effect of a laser beam perpendicular to the long axis of 
the condensate. In our model the height of the Gaussian barrier is periodic and performs 
small oscillations around its mean value. 

The Bose condensate is confined by a harmonic potential in the transverse direction while 
the external potential in the horizontal axial direction is given by: 

V(z, t) = V^t) e- z2 ' a2 + -muj 2 z z 2 , (4) 

2 

where again 

2tt 

V 1 (t)=V a + V b sin(—t) . (5) 

T 

We perform the numerical integration of NPSE with the potential (4) using the z-harmonic 
oscillator units with u z = cjj_/10 and cu± = 2n kHz. The shape of this double- well external 
potential is shown in Figure 6. Such a configuration can be experimentally obtained by 
using a procedure similar to that described in the previous section (see also [15]). 

First we calculate the "ground-state" of the system fully confined in the right well by 
choosing V\ = 20 and V 2 = 0. Then we reduce the energy barrier V 1 studying the time 
evolution and periodic tunneling of the condensate through the Gaussian barrier. In Figure 
7 we plot the fraction of 23 Na atoms in the left well as a function of time t, choosing 
V\ = 6.5 and V 2 = 0. The behavior of Phif) depends on the inter-atomic energy of the 
condensate. As shown in Figure 7, for small values of the inter-atomic strength Na s /a z 
the condensed cloud oscillates between the two wells, while for larger values of Na s /a z 
the condensate remains self-trapped. Thus, for large values of the inter-atomic interaction, 
only a small part of atoms oscillates between the two wells while the main part of the 
atomic condensate remains confined in the right well. This phenomenon, called macroscopic 
quantum self-trapping (MQST), has been recently predicted by Smerzi et al. [1]. 
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Smerzi et al. [1] have found that the time-dependent behavior of the condensate in the 
tunneling energy range can be described by the two-mode equations 

C = -v/l-C 2 sin0, <j) = A( + —j==^ cos , (6) 

where ( = (Ni — N 2 )/N is the fractional population imbalance of the condensate in the 
two wells, = 0i — 02 is the relative phase (which can be initially zero), and A = 2Ej/ E T 
with Ej the inter-atomic energy and Et the tunneling energy, i.e. the kinetic+potential 
energy splitting between the ground state and the quasi-degenerate odd first excited state 
of the GP equation. Note that E T = hu , where u is the oscillation frequency of the 
Bose condensate between the two wells when the inter-atomic interaction is zero (Ei = 0). 
For a fixed A (A > 2), there exists a critical ( c = 2V A — 1/A such that for < ( « ( c 
there are Josephson-like oscillations of the condensate with period Tj = r / y/l + A, where 
r = 2-k/uoq. But for ( c < ( < 1 there is macroscopic quantum self-trapping (MQST) of the 
condensate: even if the populations of the two wells are initially set in an asymmetric state 
(C(0) 7^ 0) they maintain, on the average, the original population imbalance with a very 
small periodic transfer of particles through the barrier with period Tst = r o/ \[%AHq — 1)) 
where H = A({0) 2 /2 - yl - C(O) 2 cos0(O) [1]. Our NPSE numerical results of Figure 7 
confirm the two-mode approximation (Eq. (6)) prediction: the fully imbalanced (C(0) = 1) 
Bose-Einstein condensate approaches a self-trapped regime by increasing Na a /a z , i.e. the 
inter-atomic energy. 

It is interesting to see what happens if we include a periodic oscillation of the barrier 
(V 2 7^ 0). As shown in Figure 8, where we plot the condensate fraction -Pl(£) in the left- 
handed well obtained by numerically integrating the NPSE, the presence of an oscillating 
barrier of period r can break or simply reduce the MQST of period t$t- In particular, in 
Figure 8 it is shown that if r = r ST j2 (parametric resonance condition [14]) then the MQST 
is broken; instead, if r = 2t S t the MQST is not broken but there is a larger periodic transfer 
of particles through the barrier. 

This phenomenon can be investigated by using the two- mode approximation of Eq. (4). 
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In this case, the presence of a oscillating barrier can be modelled by a time- dependent 
A: e.g. A(t) = A (l + e sin (2-Kt/r)). In Figure 9 we plot the population imbalance ((t) 
and its phase-space portrait (((t),((t)) with A = 25, C(0) = 0.6, 0(0) = 0. Also in this 
case the MQST is strongly affected by the parametric resonance between the period r of 
MQST oscillation and the period T of oscillation of A (i.e. the energy barrier). As shown 
in Figure 10, at the parametric resonance condition r = tst/2 and with a sufficiently large 
perturbation (e = 0.2) the system escapes from the self-trapping configuration. We have 
verified that for this system the general parametric resonance condition r = (n/2)rsT [14] 
works quite well at small n (even and odd). 

Recently the parametric excitation of cold trapped atoms in far-off-resonance optical lat- 
tices has been experimentally obtained by modulating the potential depth [16,17]. In these 
experiments the parametric resonance results in heating or losses for the trapped atoms and 
it can be studied to determine the spring constant of the periodic optical potential [17]. An- 
other kind of parametric resonance, namely the parametric resonance between the collective 
oscillations of a trapped Bose-Einstein condensate and the oscillations of the confining har- 
monic trapping potential has been demonstrated with numerical simulations by Kevrekidis, 
Bishop and Rasmussen [18]. They have shown that a weak harmonic modulation of the 
confining potential can cause an anomalously large amplitude in the collective oscillations 
of the Bose condensate. Actually, a similar effect can be obtained without parametric reso- 
nance by imposing particular values of anisotropy to the trapping potential: for such values 
different collective modes of the condensate are in resonance [19]. 

CONCLUSIONS 

We have studied the macroscopic quantum tunneling of a cigar-shaped Bose-Einstein con- 
densate confined by a harmonic potential in the transverse direction and tunneling through 
an energy barrier in the axial direction. The time-evolution of the macroscopic wave function 
of the condensate is obtained by using an effective ID nonpolynomial nonlinear Schrodinger 
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equation. 

First we have investigated a Bose condensate confined in the vertical axial direction by 
two Gaussian barriers that model a potential well. By periodically changing the height 
of the lower-lying Gaussian barrier it is possible to generate periodic waves of coherent 
matter: only when the height of the energy barrier is sufficiently small a consistent fraction 
of condensed atoms can tunnel. We show that it is possible to control the period of emission 
of the matter waves and the tunneling probability. We have also found a resonance-induced 
coherent emission of atoms: the emission probability is strongly enhanced if the period of 
oscillation of the height of the energy barrier is in parametric resonance with period of 
oscillation of the center of mass of the condensate inside the potential well. 

Then we have studied the periodic tunneling and the quantum self-trapping of a Bose- 
Einstein condensate in a double-well potential with an oscillating energy barrier. We have 
used our effective ID nonpolynomial Schrodinger equation and also the two-mode classical- 
like equations. We have found that the macroscopic quantum self-trapping of a Bose con- 
densate can be reduced and also suppressed by changing the frequency of the oscillating 
barrier. In particular, the system escapes from the self-trapping configuration if the the 
period of oscillation of the double-well energy barrier and the period of MQST oscillations 
of the condensate satisfies the parametric resonance condition. 

The parametric driving of Bose-Einstein condensates can be obtained by current ex- 
periments using optical dipole forces with far-detuned laser beams [20]. By varying the 
intensity of the laser beams one can control the height of the energy barrier that confines 
the condensed sample. The experimental investigation of the physical configurations we have 
considered in this paper may contribute to the realization of novel phenomena: controlled 
periodic emission in atom lasers and breaking of macroscopic quantum self-trapping induced 
by parametric resonance. 
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FIG. 1. Vertical axial external potential V(z) of Eq. (2). Scaled parameters of the potential: 
V\ = 300, V2 = 900, z\ = 17, Z2 = 24, a = 1.6. Length z in units a z = (h/muj z ), where uj z = ojj_/10 
with lo± = 2ir kHz. Energy in units Tiu z . 
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FIG. 2. Axial density profile of Bose condensed 23 Na atoms tunneling through the potential 
well (Eq. (5)), as obtained by solving NPSE. Time: t = 4.8. Number of 23 Na atoms: N = 10 4 . 
Scattering length: a s = 30 A. Potential well parameters: Z\ = 17, Z2 = 24 and V2 = 900. Three 
values of V\ correspond to three values of the effective chemical potential fj, e ff and of the period 
t osc of small oscillations around the minimum of the potential well: fi e ff = 53.91 and t osc = 1.21 
(top), n e ff = 54.23 and t osc = 1.17 (middle), fi e ff = 54.55 and t osc = 0.90 (bottom). Length z in 
units a z = (h/mu z ), where uj z = a;j_/10 with lo± = 2ir kHz; energy in units hu z and time in units 
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FIG. 3. Axial density profile of Bose condensed 23 Na atoms tunneling through the potential 
well (Eq. (2)) with oscillating barrier: Vi(t) = V a + V b sm(^-t), V a = 300 and V b = 200. Results 
obtained by solving NPSE. Period of oscillation: r = 1.5. Effective chemical potential of the initial 
condensate: /it e // = 57.43. Number of 23 Na atoms: N = 10 4 . Scattering length: a s = 30 A. Note 
that P T = 0.26 with Vi(t) = V a - V b = 100. Units as in Fig. 2. 
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FIG. 4. Tunneling ratio Pt/Pt at t = 5 as a function of the period r of oscillation of the 
lower-lying Gaussian barrier (Eq. (3)), as obtained by solving NPSE. V\(t) = V a + Vf,sin(^ L t), 
V a = 300 and U 6 = 200 (Eq. (6)). The period t osc of small oscillations around the minimum 
of the potential well with V a = 300 and Vj, = is t osc = 0.84. is the tunneling probability 
with Vi = V a — Vb = 100. Effective chemical potential of the initial condensate: /x e jj = 27.15 for 
N = 10 3 , fj, e ff = 57.43 for N = 10 4 . N is the number of 23 Na atoms. Scattering length: a s = 30 
A. Inset: tunneling probability Pr at t=5 as a function of r. Units as in Fig. 2. 
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FIG. 5. Axial coordinate z(t) (left) and phase-space portrait (z,p z ) with p z = z (right) of a 
classical particle under the action of the potential V(z) given by Eq. (2) with its parameters as in 
Fig. 4. From top to bottom: (a) V b = and r = anything; (b) V b = 200, r = 0.5; (c) V b = 200 and 
r = 0.75; (d) V b = 200 and r = 3. Units as in Fig. 2. 
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FIG. 6. Horizontal axial double-well potential V(z) given by Eq. (4). Scaled parameters of the 
potential: V\ = 5, Vi = 0, a = 1. Length z in units a z = (Ti/moj z ) and energy in units hu z . 
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FIG. 7. Fraction Pl(£) of Bose condensed 23 Na atoms in the left well obtained by solving NPSE 
where the external potential is given by Eq. (3) with V\ = 6.5 and Vi = 0. From top to bottom and 
from left to right: (a) Na s /a z = 0; (b) Na s /a z = 0.014; (c) Na s /a z = 0.07; (d) Na s /a z = 0.14. 
Length z in units a z = (h/mcu z ), where u> z = u±/10 with u>± = 2tt kHz. Energy in units Hoj z and 
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FIG. 8. Fraction of Bose condensed 23 Na atoms in the left well obtained by solving NPSE 

where the external potential is given by Eq. (3) with V\ = 6.5. From top to bottom: (a) V2 = 
and r = anything; (b) V2 = 2 and r = tstI% (c) V2 = 2 and r = 2tst; where tst is the period of 
unperturbed MQST oscillation. Units as in Fig. 7. 
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FIG. 9. Fractional population imbalance ((t) of the condensate in the two wells (left) and its 
phase-space portrait (C,P() with P( = ( (right), as obtained by solving Eq. (6). Initial conditions: 
C(0) = 0.6, 0(0) = 0. A = A (l + esin(27rt/r)) where A = 25. From top to bottom: (a) e = 
and r = anything; (b) e = 0.1 and r = tst/2] (c) e = 0.1 and r = 2tst', where tst is the period 
of unperturbed MQST oscillations. 
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FIG. 10. Fractional population imbalance ((t) of the condensate in the two wells (left) and its 
phase-space portrait (C,P() with P( = ( (right), as obtained by solving Eq. (6). Initial conditions: 
C(0) = 0.6, 0(0) = 0. A = Ao(l + esin(27rt/r)) where A = 25 and e = 0.2. From top to bottom: 
(a) r = t^t/2; (b) r = 2tst', where tst the period of unperturbed MQST oscillations. 
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